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Abstract 

The basic physics of bulk magnetic superconductors (MS) related to the problem of 
the coexistence of singlet superconductivity (SC) and magnetic order is reviewed. 
The interplay between exchange (EX) and electromagnetic (EM) interaction is dis- 
cussed and argued that the singlet SC and uniform ferromagnetic (F) order practi- 
cally never coexist. In case of their mutual coexistence the F order is modified into a 
domain-like or spiral structure depending on magnetic anisotropy. It turns out that 
this situation is realized in several superconductors such as ErRhiB^, HoMoqSs, 
HoMcqScs with electronic and in Auln2 with nuclear magnetic order. The later 
problem is also discussed here. 

The coexistence of SC and antiferromagnetism is more favorable than with the 
modified F order. Very interesting physics is realized in systems with SC and weak- 
ferromagnetism which results in an very reach phase diagram. 

The properties of magnetic superconductors in magnetic field are very peculiar, 
especially near the (ferro)magnetic transition temperature where the upper critical 
field becomes smaller than the thermodynamical critical field. 

The extremely interesting physics of Josephson junctions based on MS with spiral 
magnetic order is also discussed. The existence of the triplet pairing amplitude 
(-F|l) in MS with rotating magnetization (the effect recently rediscovered in SFS 
junctions) gives rise to the so called 7r-contact. Furthermore, the interplay of the 
superconducting and magnetic phase in such a contact renders possibilities for a 
new type of coupled Josephson-qubits. 
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1 Introduction 



The physics of magnetic superconductors is very interesting subject due to 
the pronounced competition of magnetic order and singlet superconductivity 
in bulk materials. The question of their coexistence was first rised theoretically 
in the pioneering work by V. L. Ginzburg [1] in 1956, where only the electro- 
magnetic (EM) interaction between magnetic moments and superconductiv- 
ity was considered. However, the breakthrough in the physics of MS came 
with experiments after the discovery of ternary rare earth (RE) compounds 
such as borides {RE)TiB4^ with transition elements T = Rh, Ir, chalcogenides 
{RE)Mo6X8 (X = S,Se), sihcides {RE)2T2Sk and stannides {RE)T^Sn^ [2]. 
In most of them type-II superconductivity is realized and in all of them are 
the localized RE magnetic moments regularly distributed in the crystal lattice. 
The basic crystallographic structure, for instance in RERh4B4, contains local- 
ized moments (LMs) which are rather far a way from the Rh and B blocks 
which deliver conduction electrons. Due to this spatial separation the con- 
duction electrons jump rarely on magnetic ions making the direct exchange 
interaction (EX) Jsf{^ 10'^ K) much smaller than in transition metallic mag- 
nets. In these systems the 4f rare-earth shells are responsible for localized 
moments in which the f-level lies much below the Fermi energy, Ef « Ep- A 
number of compounds belonging to the above families have shown coexistence 
of superconductivity with the antiferromagnetic (AF) order - antiferromag- 
netic superconductors (AFS), such as {RE)Rh4^B4^ (RE=Dy, Sm,...), and in 
most of them the Neel (AF) transition temperature TV is smaller than the 
superconducting one Tc, i.e. Tn « Tc- 

However, a lot of research, both experimental and theoretical, starting from 
the late seventies were devoted to MS systems in which ferromagnetic (F) 
order and singlet SC compete due to their strong antagonistic characters - we 
call these systems ferromagnetic superconductors {FS) and they are the main 
subject of this review. It turned out that the modified F and SC, order can 
under certain conditions coexist since the F order is transformed (in the pres- 
ence of superconductivity) into a spiral or domain-like structure - depending 
on the type and strength of magnetic anisotropy in the system [3] , [4] . In the 
RE ternary compounds this competition is rather strong and therefore these 
two orderings coexist in a limited temperature interval Tc2 < T < Tm (the 
reentrant behavior), for instance in ErRh4B4 and HoMogSs- The coexistence 
region in ErRh4B4 is narrow with pa 0.8 K, Tc2 ~ 0.7 K and = 8.7 K, 
while for HoMoeSg it is even narrower with 0.74 K, Tc2 ~ 0.7 K and 

= 1.8 X - see Refs. [2], [3], [4]. In HoMoeScg where = 5.5 K, ^ 0.8 
the exchange interaction is weaker and the coexistence persists down to T = 
K. 

A new and very interesting research field in the physics of ferromagnetic su- 
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perconductors was opened in 1997 by Pobcll's group in Bayreuth [5], which 
discovered the coexistence of superconductivity and nuclear magnetic order 
in the type-I superconductor Auln2 with = 0, 207 K and = 35 /iK. 
Although there is a tendency to the nuclear ferromagnetic order, supercon- 
ducting electrons enforce the appearance of a spiral or domain-like nuclear 
magnetic ordering in the SC state below [6] . 

It turns out that not only bulk properties of FS are exotic, but also Josephson 
junctions made of bulk MS with spiral ordering show potentially fascinating 
properties, such as 7r-contact [7], combination of a magnetic analog of the 
Josephson effect for spin current with the ordinary Josephson effect for charge 
current [8]. 

In the following we shall discuss mainly the microscopic and macroscopic the- 
ory of ferromagnetic and antiferromagnetic superconductors which takes into 
account the most relevant interactions between localized moments and con- 
duction electrons - the exchange (EX) and electromagnetic (EM) interaction. 
This theory was elaborated by Buzdin, Bulaevskii, Panyukov and present au- 
thor - see [3] and references therein, and successfully applied to a number of 
systems. Due to the lack of space we shall discuss effects of magnetic field 
on magnetic superconductors briefly - for this subject we refer the reader to 
Ref.[3]. 

We would hke to point out here that in the last several years there was a huge 
activity in studying of hybrid heterogeneous magnetic superconductors such as 
S-F multilayers and S-F-S Josephson junctions. This field is not only of impor- 
tance for the fundamental solid state physics but it is of enormous interest for 
applications in spintronics and quantum computing, especially after the exper- 
imental confirmation of the remarkable prediction of the tt- Josephson contact 
by Alexander Buzdin and coworkers [9], [10]. This exciting field will be covered 
elsewhere in this issue, as well as the physics of other magnetic superconduc- 
tors - heavy fermions, borocarbides {RE)Ni2B2C, cuprates RuSr2GdCu208, 
ferromagnets with triplet SC such as UGe2. 



2 Competition between SC and F order in FS 

Here we shall be limited to those magnetic superconductors where the mag- 
netic ordering of the localized 4f moments (LM) is due to the indirect exchange 
interaction (RKKY) going via the conduction electrons. The characteristic ex- 
change energy is 0,.^. ~ N{0)hlj where iV(0) is the density of states at the Fermi 
level (per LM) and ho{= {g — l)nJsf{0){Jz)) is the maximal exchange field act- 
ing on conduction electrons. Here, g is the Lande factor, n is the density of lo- 
calized magnetic moments (LMs), Js/(0) is the direct exchange energy between 
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conduction electrons and LMs, (J^) is the averaged total angular moment of 
the LM. Note that the exchange field acting on electrons is hex{r) = hoS{r), 
where S{r){= {Jz)/J) is the localized spin normalized to one. Let us mention 
in advance that in the RE ternary compounds the exchange field Hq is still 
rather large, i.e ho ~ 10^ K and ho ^ Aq < 10 K. We shall see below that 
in spite of the fact that ho is larger than the Clogston paramagnetic field hp, 
i.e. ho ^ hp ^ 0.7Ao, there is a coexistence of SC and modified ferromagnetic 
order. In the presence of magnetic ordering characterized by the magnetiza- 
tion M(r) there is electromagnetic interaction between localized moments and 
(super)conducting electrons, since M(r) = n/iS(r) creates the dipolar mag- 
netic field B(r) = rot A{r) which on the other side induce screening current 
js of conduction electrons (the Meissner effect). The Fourier transformed is 
related to A by the kernel i^s(q), i.e. js(q) = — -ft's(q)A(q). Having in mind 
magnetic superconductors based on RE ternary compounds we shall discuss 
the physics in the mean-field approximation for SC and magnetic subsystem 
- quite appropriate approach, i.e. we put S(r) — > S(r) =< S(r) > 

The total Hamiltonian of the system is given by 

H = Hbcs + Hcf + H^^p + E[-B(ri)y//BJi + HcF{Ji)] (1) 

i 

Hbcs = J ciV{V^t(r)[e(p-^A)] - /.]V;(r) + V^t(r)v;,(r)V^(r) 

+ lA{r)i;^{r)iaJ\r) - lA*{r)i^{r)iaJ{r) + [A^}. (2) 

^ 9epi 



Here, e(p— ^A) is the band energy of electrons in magnetic field, A(r) is the 
singlet superconducting order parameter, a = {a^, ay, a^} are Pauli matrices, 
while the exchange field acting on electrons is given by 



VEx{r)-- 



\ 



(3) 



We go slightly in advance by informing the reader, that in the SC state the 
ferromagnetic order is modified into a spiral or domain-like structure with 
the wave vector Q, depending on magnetic anisotropy described by Hcf- If 
the magnetic anisotropy is small (or easy plane) than the spiral structure is 
realized with h-^{r) = he^^^ and h^{v) = 0, while in the opposite case with 
an easy axis anisotropy one has h^{r) = and h^{r) = /i^(r + L), L = 2tt/Q. 
The effect of nonmagnetic impurities is described by Himp whose effect is 
characterized by the mean-free path I and time r. 
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2.1 Sinus magnetic structure due to SC for T < 

In the RE ternary magnetic superconductors in which the singlet SC and 
ferromagnetic order compete, the superconducting critical temperature, Tc, is 
much higher than the magnetic one, i.e. << T^. Before discussing the com- 
plete phase diagram we shall study the coexistence problem at temperatures 
near T^, i.e. T < T^, where the magnetic order parameter is small S << 1. 
In case when the easy-axis magnetic anisotropy D is sufficiently large then 
the sinus structure S(r) ^ S{T) sin Qr appears below Tm (for small D a spiral 
order is favored - see 2.3). In that case h^xi^) = ho \ S(r) |<< ho, A and the 
free-energy can be calculated by the perturbation theory 

F{S(r), A(r), A(r)} = Fm{S} + Fs{A} + F,„,{S, A, A}, (4) 



where Fm and Fs are the magnetic and SC functional without mutual inter- 
action,respectively. 

FM{S(r)} = - T^o) + ea^)] I S(q) \' -D \ S^) ?) 

q ^ 

/■^3 (B-47rM)2 

where S{q) is the Fourier transform of S{r). The last term in Eq. (5) is the 
magnetic energy for a given magnetization M(r) = n/iS(r) and the magnetic 
induction B = 47rM. The characteristic energy for the EM interaction is given 
by 9em = (B'^/Snn) = 27rn/x^ which is ~ 1 in the RE ternary compounds. 

Since T^KKT^ and /i < A the SC free-energy density (F5 = / (frPs) 

F5{A(r)} = -lAr(0)AMn^ (6) 

is minimized for A Aq and we omit it from the analysis near T^. The part 
Fint describes the EX and EM interaction between SC and magnetic order 
(note j,(q) = -ii',(q)A(q)) 

Fim = Y.{\Ks{ci) I A(q) ^ +^e.^^^^^^^7^ | S(q) 1^. (7) 
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After minimizing F{S(r), A(r), A(r)} with respect to A(r) one gets F{S(r), A(r)} 
in the following form (see more below and in [3]) 

F{S(r),A(r)} = n^{^[(T-7;„o) + ^aV)] I S(q) ^ -D. \ SM) ? 

q ^ 

^ Oex Xn(q) -X.(q) I . ,2 , e.m ^'KK.^q) I S(q) p +(qS(q))(qS(-q)) ^^^^ 
+ T Xn(0) I ^^^^ I +^ ,- + AnKM) ^^'^ 



Here, a is of the order of lattice constant (magnetic stiffness) and the bare 
critical temperature T^q and 9 take in a subtle way (note T^o 7^ ^ - see details 
in [3]) into account the indirect EX and direct dipole-dipole (EM) interaction 
between LMs - see [3]. Xn(q) and Xs(q) are electronic susceptibilities in the 
normal and SC state, respectively. 0,,,,, = 27ifi'^ characterizes the EM effects 
in the dipole-dipole interaction between LMs, while Kg{q) is the EM kernel 
which describes the screening effect of the dipole-dipole interaction by the- 
superconducting electrons. Dz{> 0) is the magnetic anisotropy which orients 
spins along the z-axis. 

Due to the singlet SC pairing Xs(q) is reduced significantly at small wave vec- 
tors q < where is the SC coherence length. In the clean limit {I — > oo) 
and at T = one has Xs(0) — which means that the ferromagnetic order can 
not coexist with singlet superconductivity. In Fig.l we show Xs,n(q) schemati- 
cally for the cases when the ferromagnetic (la) or antiferromagnetic order (lb) 
is realized in the normal state. It is seen that a singlet superconductor behaves 
as a normal metal at large momenta, i.e. Xs{q ~ kp) is weakly affected by SC. 
Therefore AF competes with SC much less than the ferromagnetic order does. 

We stress that at finite temperature Xs(0) ^ but exponentially small in 
singlet s-wave SC, while in d-wave SC one has Xs(0) ~ x„(0)(T'/Ao). In the 
presence of the spin-orbit (SO) scattering Xs(0) is also finite. The general 
expression for Xs(q) is calculated in [11] 



where — cu/A and 



Pi^, q) = I (10) 

2 arctan{^ Jl + ul + l/2r_} 



Here, ^ = -{4:/3)TgJ, tJ=t ^ + Tg^, Igo = vpTso anda;„ = 7rT(2n-Hl). 
Later, we shall discuss some effects of the SO interaction on the coexistence 
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Figure 1. Schematic spin susceptibility in the SC and normal state Xn,s{fl)' (a) for 
the ferromagnetic order in the normal state peak at Q = 0; (b) for the antiferro- 
magnetic order - peak at Qo- 



phase. The effect of the exchange scattering is similar, i.e. Xs(0) is finite for 
finite r^. 

Since in the following we study the competition between SC and ferromag- 
netism at low temperatures it is sufficient to give a general expression for 
Ks{q) in the clean limit 

mcqvp J X Jl + (3|Ei:)2 



The expression for finite / is cumbersome and is omitted here. Some limiting 
cases of Ks{q) which are relevant for real magnetic superconductors will be 
studied below. 

By minimizing F{S(r), A(r)} w.r.t. the wave vector q one obtains the equilib- 
rium magnetic structure which depends on microscopic parameters a, .^O) AL(the 
London penetration depth), 9ex,0em- From the above equation we conclude 
that due to the EM interaction the magnetic structure is transverse, q • S(q) ~ 
0. 

Let us analyze ^^^id Ks{q) in the interesting range of parameters. In the 

clean limit and for q^g ^ 1 one has 

Ks{q) = Xn(q) - X.(q) = Xn(0)(l - ^^), (12) 



7 



while for g^o ^ 1 it holds 

^^^^^ = Zx^;^"^^) - = ^1^) 



Based on Eqs. (12-13) and the expression for the free-energy F in Eq.(8) we 
obtain that just below the transition temperature T,„ = Tmo(l — 3(f %^|^)^^^) 
a transverse (Q^XS'^) sinus structure Szijc) ^ 5'sin(Q5r) appears. In case 
when -C aX^ the wave vector Qg is determined by the EX interaction and 
for Oejeem > (a/Co)' it is given by [12], [13] 



For 9ex/^em ^ (t^/'Co)^ the EM interaction prevails 

Q. ^ (^)'/^. (15) 

In the opposite limit ^ aX^ the EX interaction dominates for 9ex/0em ^ 
{a'^Co/Xlf/^ which gives again [12], [13] 

Q. = (J^)"^. (16) 

For 9ex/0em ^ ("^^Co/Al)^/^ the EM interaction dominates which gives 



From these expressions it is seen that in realistic cases Qg is determined by 
the EX interaction - it does not depend on the EM parameter A^, while the 
EM interaction (with dependence of Q) is dominant only for extremely 
small EX interaction {O^x ^ ^em(a/Co)^ or (^ex ^ ^em(ct^Co/'^i)^^^); i-C- for 
{dex/dem) ^ 10^"^ — 10^^. Howcver, in typical ferromagnetic superconductors, 
such as ErRh4B4, HoMoeSs, HoMoeSeg, Auln2, the EX interaction dominates 
since 9ex > 0.1 9em and a <^ < Xl- 

In reality nonmagnetic impurities are always present and one should know 
Xs(q, and Kgifii I) as a function of the mean-free path /. The corresponding 
calculations show that if {l^ /a^^^X^) <^ 1 one has for 9^x1 Oem ^ «^Co/^^ [12], 
[13] 
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and 



for {a%/l') » 0e,/eem » (^7^1), aH/Cl 

In case when {9 ex /O em) ^ (a^Co/^^) or {9 ex /O em) -C ^^/''^i the EM interaction 
dominates 

Qs - i^f"- (20) 



Let us stress some interesting properties of ferromagnetic superconductors: 
(i) the ferromagnetic critical temperature is strongly reduced in the presence 
of SC due to the formation of Cooper pairs in the SC state, i.e. one has 
Tp = Tmo(l — issd^sm.^ <^ x^- In fact this result is more general and holds also 
for the coexistence of SC and itinerant ferromagnetism (F) - singlet SC and 
ferromagnetism do not coexist. In that sense a number of recent papers which 
claim that the itinerant F and SC coexist should be completely abandoned 
[14]. However, in some itinerant ferromagnets such as I9C07 (with Tp = 4.5 
K) the microscopic parameters favor spiral or domain magnetic structure in 
the SC state with Tc = 2.5 K as it was proposed in [15]; (ii) in isotropic 
magnetic systems and near the critical temperature the inverse scattering 
time due to magnetic fluctuations can diverge and thus destroy SC. How- 
ever, this divergence is suppressed in the real RE ternary compounds due to 
the long-range dipolc-diploe interaction. The interaction of SC with magnetic 
fluctuations is described by the free-energy contribution 

\Xn(q)-X.(q) 

Fscji - -Y :^{Sz,,Sz,.,) , (21) 



where 

T + a-'q^ + {eem/e)ql/q 



{Sz,qSz,-q) ^ , „2„2^(a //)N^2/^2' (^2) 



with T = (T— Tmo) / 0. Due to the large dipole-dipole temperature with 6em ~ ^ 
these fluctuations looks like four-dimensional, thus giving rather small value 
for the inverse scattering time ~ ^ ^ Tc; (iii) the relative strength of the 
EX and EM interaction is controlled by the parameter 

_ ^Int _ ^em ^ /oo\ 

Ff^f^' Q.xQ^\l ^ ^ 
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Figure 2. The striped domain magnetic structure S{x) = Sz{x)ez with the period 
Ld = 27r/QD5; Q is along the x-axis. 

In the RE ternary compounds the case r <^ 1 is always realized, due to the 
large value of Q^\\ ^ 1. Therefore, practically in all RE ternary compounds 
the EX interaction dominates in the formation of magnetic structure, while 
the EM interaction makes it transversal - see exception in weak ferromagnets 
below; (iv) In spite of the fact that in the RE ternary compounds the ferro- 
magnetism is stronger phenomenon than SC - the gain in the ferromagnetic 
energy (per LM and at T=0 K) Em ~ N{{))h'^ is larger than the gain in the 
SC condensation energy ^ A^(0)Ag since /io(~ 10^ K) > Ao(< 10 K), the 
ferromagnetic order is more "generous" and varies spatially in the SC state. 
The reason lies in the fact that the magnetic stiffness (~ a) is much smaller 
than the superconducting stiffness (~ ^o)- 



2.2 Domain magnetic structure due to SC 

By lowering T the higher term ~ S'^{r) makes the change of the moduo of S(r) 
unfavorable and the sinus-structure is transformed, as it will be shown below, 
into the striped domain structure (DS) - see Fig.2. At the same time since the 
exchange field grows h^x = hoS(T) but for hexiT) < A the mutual interaction 
of magnetism and SC can be treated by the perturbation theory. In such a case 
the free-energy density is completed by the density of the domain wall energy 
QEw /t^, where Ew is the wall-energy per unit surface. In case of sufficiently 
large magnetic anisotropy > 6 rotation of the moments in the wall is 
unfavorable and the linear domain wall is favored, i.e. S'^(x) = Sth{x/lw), 
Sx = Sy = 0, where Iw = aj ^ is the domain-wall thickness [3]. The domain 
wall energy per unit surface is given by 
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Ew = {AV2l?,)neS'^aT^''^ = nOS'^a (24) 
where r = {T — Tjno)/9) - see [3]. The free-energy density Fds in the DS phase 



IS 



Fns = neiKs' + ^-S'] + ^E^ + ^^e.^^- (25) 

Z 4 TT ZTT C^^o 



Note, that if the anisotropy energy is small, i.e. r > 20^/6 the rotating domain 
wall is realized with the wall thickness l^r ^ a{QjD^^I'^ and the wall energy 
E-i^ ~ S'^a{ODzY^'^ . In the case of an extremely small anisotropy (D^/^) < 
(a/^o) ~ 10~^ — 10~^ then the spiral structure is realized. Minimizing Fds 
in Eq.(25) with respect to Q one obtains the equilibrium wave vector of the 
striped DS phase 

Qds ^ (26) 

We conclude that the period of the DS structure is larger than for the sinus 

phase. It is worth of mentioning that: (i) the here obtained striped domain 
structure is due to SC and it is property of the bulk, (ii) the problem of the DS 
phase in SC is mathematically similar to the problem of domain structure in a 
normal ferromagnetic plate with the magnetization perpendicular to the plat 
plane, where the role of Fj„( is played by the magnetic energy dissipated out 
of plate. Generally, the domain structure is realized when the wall thickness 
o/a/t is much smaller then the striped domain thickness tt/Q implying that 
r » (a/Co)'/' ~ 10-'. 

At lower temperatures when hex{T) > A the nonperturbative problem is stud- 
ied in the presence of impurities by the quasiclassical ELO equations (see Ap- 
pendix 7.2) since the period Ld oi the domain structure is much larger than 
a, i.e. Ld 3> a. This equations are solved for the domain structure with 

^ fe=0 z/e -h 1 q 



By assuming that Q is along the x-axis the solution for the Green's function 
are searched in the form 

/(v,x) = /o(v) + ^/,(v)e^^ (28) 

k 



and analogously for g((v,a;), where k = (2m + 1)Q and | /*; |<C| /o |, | S'fc |^| 
go \. The calculations were done in [12], [3] and here we present only the final 
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result for the free-energy in the dirty limit {I ^ ^o)- It turns out, that in that 
case the interaction of the magnetic domain structure with SC is similar to 
the case of magnetic impurities with the inverse scattering time and with 
TmA > 1, i.e. Fds is given by 



FDs^n9[^TS^ + ^S^]+QEw 



l^(0)A^ln(f) + ^(0)|^(l-^). (29) 



is given by {h^^q = hoS^^q) 



where 



L 2y arctany 

^ 7r(y — arctany) 



and 



L2(2/) = -[(l + 4)arctany--]. (32) 
TT y^ y 



The magnetic induction is given by 

(33) 



A7mfi[q'^Sg - q(qSj] 



g2 + X,(5)(l-4/37rT„A)' 



where the Kernel Kg(q) in the dirty limit has the form 



Based on the free-energy in Eq.(29) we can study the coexistence problem 
in the whole temperature regions and for various Ql - see more in [12]. We 
summarize the main results: (i) at T = the sinusoidal magnetic order 
appears with the wave vector Qg ~ (1/a^^o)^^^ ! (H) by lowering temperature 
the domain structure appears with Q^s ~ (l/'^Co)^''^ which persists up to the 
temperature of the first order phase transition Tc2 where the DS phase passes 
into the normal ferromagnetic state. At Tc2 one has 

Fds{Sds{T,2), A(T,2), Qds,c2} = Ffn{Sf{T,2), 0, 0}, (35) 
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where Qds,c2 ~ 1.8(a(T,2)eo)-'/' ~ «o)-'/^ A(T,2) = O.SSAq and {S!,/Qc2) 
0.07(Aot;F//io); if Sds{Tc2) > 1 then the DS is stable up to T = 
X - this situation is reahzed in systems with small EX interaction (which 
still dominates over EM), i.e. for 9ex < Ogx ~ (^c/^o); i'^''^) dirty SC with 
{hrY <S 1 there is a gap in the quasiparticle spectrum for £^ < A in the whole 
range of the existence of the domain phase. The calculations in clean SC show 
[4], [3] that when h(T) ^ A the spectrum is gapless. For instance in the DS 
phase one has for £^ <^ A 



N{E) _ Tih E 4A 
N{0) ~ ^A 

while in the case of the spiral order 

N{E) _ nh E 
N{0) ~ 2vfQA' 



(36) 



(37) 



(v) The spin-orbit interaction decreases the value of x„(0) — XsW (and small 
q) and it is detrimental for the DS phase. However the analysis in [16] shows 
that the S-0 scattering destroys the peak in Xs(q) only in very dirty systems 
when / ~ a. 

We would like to point out, that there were a lot of studies of ferromagnetic su- 
perconductors based on the phenomenological theory which takes into account 
the EM interaction only [17], [18]. However, this interesting phenomenology 
is inadequate in describing real materials, such as the above numbered RE 
ternary compounds where the EX interaction prevails in the formation of the 
oscillatory structure (with Q3> XJ,^ ) in the SC state. We stress that in 
the above theory the EM interaction, as well as the EX one, is taken into 
account on the microscopic level, thus giving much more reliable predictions 
than the phenomenological approach. 

2.2.1 Experimental situation 

In the most important ferromagnetic superconductors HoMogSg, ErRh4B4, 
HoMoeScg the range of microscopic parameters allows the coexistence of SC 
and modified ferromagnetic order. For instance in clean systems one has [3]: 
in ErRh4B4 - n ~ 10^^ cm'^, /x = 5.6 /^b, a ^ 1 i, Al(0) « 900 i, ~ 200 
i, Ao ^ 15.5 K, Ar-i(O) = 1850 K ■ spin RE, vp ~ 1.3 x 10^ cm-^s, 6^^ ~ 0.5 
ho ^ 40 K, T-^ ^ 3 K and e^m ~ 1.8 in HoMogSg - n - 4 x 10^^ 
cm"3, n = 9.1 /i^, d ^ 2.5 i, Xl{0) ~ 1200 i, ~ 1500 i, Aq ^ 3.2 K, 
N-\0) = 3600 K ■ spin RE, vp ^ 1.8 x 10^ cm-\s, O^x ^ 0.2 K, ho ^ 24 K, 
'^rn ~ 0-9 K and 6em ~ 1-3 K while in HoMoeSes a number of parameters are 
similar to HoMoeSs. 
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An oscillatory magnetic structure (either sinus or domain-like) due to SC has 
been observed at least in three compounds: (1) in HoMogSg where = 8.7 
K, T„ ^ 0.8 K, f« 0.7 K and Qds ~ 0.03 i"^; (2) in ErRh4B4 where 
= 8.7 K, ~ 0.8 K, %2 ~ 0.7 K and Qds ~ 0.06 i'^; (3) in HoMogSeg 
where = 5.5 K, = 0.5 K and Qds ~ 0.09 - 0.06 i"^ and the coexis- 
tence persists up to T = K\ All these results are in a satisfactory agreement 
with the above theory. The conclusion is that in most RE ternary compounds 
the EX interaction is responsible for the formation of the oscillatory mag- 
netic structure in the SC state, while the EM interaction makes the structure 
transverse Q • S(r) = 0. 



2.3 Domain magnetic structure in thin SC film 



In the above calculations we have assumed that the thickness L of the sample 
is very large, i.e. L 3> so that the dissipated magnetic energy (stray field) 
can be neglected. In case of thin films with L C,o the stray magnetic energy 
Est existing around the domain walls must be added to the free-energy Fds 
in Eq.(29) (or its simple version in Eq.(25)), i.e. F^ot = Fds + is given by 
[13] 

Ftot/n = {FdsM + Est = {FdsM + O-^S^em-^- (38) 



In case when the ratio r{— Fj^^^'/Fj^^'^) <^ 1 the minimization of Ftot w.r.t. 
Q gives 

Qlt-Qls + Ql. (39) 



where Qds is the wave vector of the DS phase without stray magnetic energy 
and Qf ~ ^■Q{Oem/(^ciLy^'^ is the wave vector of the domain structure in 
the normal ferromagnetic state. Prom Eq.(39) it is seen that the period of 
the DS {d — 2T^/Qtot) in thin SC film is decreased due to the stray field. It 
comes out from Eq.(38) that the transition temperature Tc2 (for the first order 
phase transition DS FN{domain) can be pushed to zero when L < = 
3^Q{6ernS^2i^ = oo))/6'ea:(l — S^2{^ = '^))^- The experiments on thin films of 
HoMoeSs show such a thickness dependence of Tc2 where Tc2{L) < Tc2(oo). 

Let us mention that even in the normal ferromagnetic state, which is real- 
ized for T < Tc2, there is possibility that SC exist in the domain walls as 
it was shown in [19], [20], [21], [3]. This situation can be realized in some 
pseudoternary compounds where /iq ^ ^o- 



14 



2-4 Coexistence of Nuclear Magnetism and Superconductivity 



In 1997 the Pobell's group from Bayreuth made an important discovery [5] 
by observing that superconductivity and nuclear magnetism coexist in Auln2 
with Tc = 0, 207 K and = 35 fj^K. At first glance this is not too surprising 
having in mind smallness of the hyperfine interaction between conduction elec- 
trons and nuclear spins. However Buzdin, Bulaevskii and the present author 
applied in 1997 [6] the theory of magnetic superconductors [3] and found a 
surprising result - the effective nuclear "exchange" field (the hyperfine contact 
interaction) is rather large hhyp ^ 1 K while Aq ~ 0.6 K, i.e. hhyp > Aq! The 
hyperfine interaction has the same (mathematical) structure as the exchange 
interaction between the 4f LMs and conduction electrons 

He-nuc - jd^^Y^ Ahyp^i^ ' R^)^|J\r)aii^|J{r) (40) 



Here, Afiyp is the hyperfine interaction and the "exchange field" is given 
by hhyp = nAhyp{ Ij), where Ij is the nuclear spin. So, the nuclear mag- 
netism in Auln2, which shows strong tendency toward ferromagnetism, com- 
petes rather strongly with SC. It was estimated from the experiment [5] that 
0,m{^ 2nnnfJ,l) ~3 l^K and 0ea,{^ N{Q)hU ^ 35 fxK, ~ 10^ A, \l ^ 10^ 
A, Z 3.6 X 10^ A {I < ^o) which means that the "EX" (contact) interaction 
is much stronger than the EM (dipole-dipole) one and the theory invented 
for RE ternary compounds is completely applicable to this case. This theory 
predict, that if the nuclear magnetic anisotropy (due to the dipole-dipole inter- 
action) is small, i.e. {D/9ex) < 10~^, the spiral magnetic structure should be 
reahzed, while in the opposite case {D/9ex) > 10~^ the striped domain struc- 
ture is formed. The experiments in magnetic field [5] give evidence that SC 
and oscillating magnetic order coexist up to T = fT, i.e. the case 6ex < ^Ix is 
realized in the type-I superconductor Auln2- Unfortunately, until now there 
were no nuclear scattering measurements on Auln2 which could resolve the 
nuclear magnetic structure below — 35 ^K. 

The study of the coexistence of SC and nuclear magnetic order is of enormous 
importance for the fundamental physics. These systems give an opportunity 

to study the coexistence problem in cases when the electronic temperature 
(Te) is different than the nuclear one (T„), i.e. Tg 7^ T„. But probably the 
most interesting problem is the coexistence of SC and nuclear magnetism in 
the case of negative nuclear temperatures {T^ <Q K). 
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3 Antiferromagnetic superconductors (AFS) 



3. 1 Coexistence of antiferromagnetism and superconductivity 



An evident experimental fact in the RE ternary compounds is that super- 
conductivity coexists with the antiferromagnetic (AF) order much easier than 
with the modified ferromagnetic order. The reason is that the effective ex- 
change field in AFS varies on the lattice constant (the AF wave vector is 
Qaf ~ and it is averaged to zero over the volume of the Cooper pair ^q- 
Thinking in terms of the electronic susceptibility one has 



Xn(0) vfQaf Ef ' ^ ' 



which means that in AFS is very small. Due to the same reason the 

EM interaction is small since SKs{Qaf) '""^ '5'^/(^i'Co)i i-C- Fj^^^^\<^ Fint^'') 
This result is confirmed in a number of RE ternary compounds in which the 
Neel temperature T/v(~ A(0)/i^) is in most cases (much) smaller than [2]. 
In these systems the magnetic scattering above Tjv is not harmful for SC since 
the inverse life time ~ T^v <S T^. 

However, there are a number of interesting properties of the AF superconduc- 
tors (AFS) such as the pair-breaking effect of nonmagnetic impurities charac- 
terized by the life-time r. In case when <^ T^. the nonmagnetic impurities 
in the AFS state are pair-breaking, like magnetic impurities with the inverse 
scattering time 

= (42) 



For hr <^ 1 one gets ~ Tjy -C Tc which means that in this case the pair- 
breaking effect of impurities is rather small [22]. Very interesting situation 
appears for systems with T/v ^ Tc. Even in such a case the exchange field 
docs not suppress Tc significantly since the theory (based on Eqs. (63-64) in 
Appendix 7.1) predicts that {5Tc/Tco) ~ {h/ Ep)(\nh/ Ep) <C 1. However, in 
the presence of nonmagnetic impurities T^ is renormalized appreciably and SC 
disappears for I < iO^oih/vpQAF) ~ Co^jv/Zi- In that respect there is one 
very interesting AFS compound Th2Mo^S4^ with TV = 19 X and T^ = 0.8 K 
where one expects that SC should disappear due to the strong magnetic scat- 
tering. However, it turns out that in this system the magnetic anisotropy, in 
conjunction with large momentum J = 9, strongly suppress this pair-breaking 
effect giving rise for SC. 
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3.2 Weak ferromagnetism in AFS 



In the case of competition of SC and the ferromagnetic order in the RE ternary 
compounds the theory predicts that in the presence of an appreciable EX 
interaction SC can coexist only with spiral and DS (or sinus) order - depending 
on the magnetic anisotropy, while other phases are excluded. It turns out 
that in AF superconductors with weak ferromagnetism (WF) - of the Moriya- 
Dyalozhinski type, the phase diagram can be much reacher containing also 
the Meissner phase {M 7^ 0, S = 0) and the spontaneous vortex state [23]. 
We discuss this problem briefly by studying the AF order with two sublattice 
where the AF order parameter is given by I = Si — S2. In systems which allow 
WF there is an additional term in the free-energy Fy^p = D[Si x S2] in the 
total free-energy. If for instance I is along the xy- plane and D is so oriented that 
it allows the appearance of the weak ferromagnetism m = S1-I-S2 (M = n//m) 
in the xy-plane then Fwf is given by 

FwF = f3n6ex{mjy + ruylx). (43) 



Since in most systems m ~ 10~^/ it immediately implies that /5 <^ 1. In 
that case and when <^ the interaction part of the total free-energy 
{F — Fjn+Fs+Fint) is given by Eq.(7), while the magnetic system is described 
by 

Fm^ I d^rneexiaif + ^{1^ + bm" + a^Vlf] 

+ J d^r[pnee:c{mjy + rUyQ + - 4jrM) ^^^^ 



By minimizing F w.r.t. A, I, m and q we get possible phases in AFS with 
WF [23]. The resulting free-energy is similar to the case of ferromagnetic 
superconductors with an effective magnetic stiffness Oe// = (ab/P) 3> a. It 
turns out that ii f3 ^ a/C,o the EX interaction dominates in the formation of 
the magnetic structure, and the sinus structure (/ ~ sin Qr and m ~ sin Qr) 
appears at T^v, while for (a/XL) < P <C a/Co the EM interaction prevails in 
the formation of the sinusoidal structure. If /3 < {a/ XL)\J'29gm/dex than the 
nonuniform structure is unfavorable and the so called Meissner state (first 
proposed by Ginzburg in 1956) appears. It is characterized by M =const and 
B = in the bulk sample due to the screening SC current on the surface 
{B = AnM exp{—z/ Xl}). By lowering the temperature | 81^2 | grow and it 
is necessary to take into account higher order terms in F. As a result one 
gets that for (3 » y^a/Co again the EX dominates and the striped DS phase 
is realized, while for Ja/XL <^ P <S v^/^o striped DS phase is realized 
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due to the EM interaction. However, by lowering the temperature the domain 
wall energy grows and it my happen that a spontaneous vortex state (with 
AnM > Hci - the lower critical field) appears for jS ^ y^a/^o and for the AF 
vector / > /, ~ (i/,i/M(0))(/?Ai/52)i/:\ 5, = a[{TN - T)/TnY/\ From the 
known RE ternary compounds the good candidate for such a behavior is the 
body centered tetragonal (b.c.t.) system ErRhiBi. 



4 Magnetic superconductors in magnetic field 

There are a number of interesting effects of the magnetic field H either in 
the coexistence phase or above the magnetic transition temperature where 
S{T) — 0. We discuss some of them briefiy - for more details see [3], [16]. 

(1) DS phase in magnetic field - In the case of bulk samples magnetic field 
penetrates only on the length Al, thus affecting surface of the sample only. 
However, in thin films the paramagnetic effect of the field is most important 
[3]. This problem was studied in the case of a thin (along the y-axis) film 
the thickness Ly < when the magnetic field is parallel to the domains, 
i.e. H = - see Fig. 2. It is found that the magnetization Sz{x) contains 
besides the odd harmonics also the zeroth-one as well as the even ones 

00 o c 

SJx) ^S5 + y —{[1- (-l)'=cos(7rM)] sm( kQx) 

+ (-1)*^ sin(7rM) cos(A;Qx)}, (45) 

with 5 = iiH/2S9ex- This change of harmonics in Sz{x) can be observed by 
magnetic neutron diffraction experiments. The result in Eq. (45) means that 
the domains with M parallel to H increase their thickness d d{l + 5) while 
those antiparallel decrease it d — > d{l — 5). In case when the zeroth component 
of the exchange field h{= hoS5) >K = A[l - (1/t^A)2/3]2/3 the DS phase is 
destroyed due to the Zeeman effect making A = 0. For h < he the parameters 
of the DS phase are renormalized, i.e. Q{H) < (5(0). In case when H = Hey 
(i.e. orthogonal to the z-axis) then all domains have the same thickness and 
there is no redistribution of intensities of neutron peaks. However, there is 
only a decrease of intensities of (2k+l)Q peaks by the factor (1 — S'jJ where 
S± = fiH/S{6ex + Dz) and is the magnetic anisotropy. 

(2) MS in magnetic fi.eld at T > Tm - The effect of the exchange field on SC 
in magnetic field is negligible for T < Tc since for Tm <^ Tc the magnetic 
susceptibility Xm is very small. However, at T near Tm there is a significant 
increase of Xm a-nd accordingly the increase of the paramagnetic effect. 
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(i) Thermodynamic critical field Hc{T) - We illustrate this effect by ana- 
lyzing the change of the thermodynamical field Hc{T) (for the transition 
N —>■ MS) in magnetic superconductors. In that case the Gibbs energy den- 
sity of the paramagnetic normal phase is equal to that of the SC phase, 
Gn{Hc) = Gsc{Hc) where 

GscW = F„(0) - ^ (46) 
Gj,{H,) = F„(0) - ^ (47) 



which gives the critical field 

H^(T) = '''f^\ (48) 



Here, {Hl^/M) = N(0)A'^/2 is the SC condensation energy and the mag- 
netic permeability is = 1 -|- 47rXm (here we neglect the conduction elec- 
tron susceptibility Xe since in MS one has Xe ^ Xm)- For T > one has 
Xm{T) ^ {9emTmo/^T^0)/{T - T^q) and 

H,{T) ~ - T^o (49) 



i.e. Hc{T) is drastically reduced near T^o- 

(ii) Upper critical field Hc2{T) - Above T^o in the presence of the external 
field He superconductivity is suppressed by the orbital effect of the field B = 
Hj(l -I- 47rXm) and by the paramagnetic effect of the exchange field h (and 
by the much smaller effect due to B). Here, Hj = Hg + H/) where H^j is the 
demagnetization field. The critical field can be calculated by the same formula 
as for usual SC - see [24], where /^s is replaced hj p,B = IJ^b + hoM/niiHi and 
the electron charge e by e = e(l -|-47rM/ifj). In the pure hmit and for T <^Tc 
one gets the modified Gruenberg-Giinther formula [25] 



where H*2{0) us the upper orbital critical field in absence of magnetic moments 
and f{a) is calculated numerically in [25]. The parameter a describes the 
relative role of the orbital and paramagnetic effect 

2H:,iO)hoXrniT) 

{1 + A7rXm{T))n^iAo ^""^ 
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In the RE ternary magnetic superconductors one has ho ^ Aq and n/i is one 
order of magnitude smaher than H*2{0) thus giving ct 3> 1 in the region where 
T <C Tc- It is known that for a > 1.8 [25] in pure superconductors it is realized 
the Larkin-Ovchinikov-Fulde-Ferrell (LOFF) phase (due to paramagnetic ef- 
fects) where the SC order parameter oscillates being also zero at some points. 
For a » 1 one has f{a) 1 and 

H,2{T)^1.5^^{T-T^o) (52) 
rloljLu 



i.e. Hc2{T) depends linearly on T — T^o near T^o, which is much faster falloff 
than of Hc{T). This leads to an very interesting effect that by the first order 
phase transition at Hc{T) the system goes into the Meissner or vortex state 
depending on the relation between He and Hd- Let us mention that the lower 
critical field Hd is very weakly affected by the exchange field due to the 
localized moments. The theory [3] predicts the following dependence of Hd 

i/ci = ^ln^, (53) 



where p takes into account screening effects due to EX and EM interaction 
p{T) = 1 ^'^'"..^ ^ (54) 

"em ~r ^ex ~r j'^^ 



Note, the theory based on the EM interaction (which assumes Q^x = 0) gives 
p(T — * Tmo) = 0, which makes the effective penetration depth zero, i.e. Ag// = 
\l\/p = 0, and the Ginzburg-Landau parameter k = (Xeff/C) 0. If this 
would be correct than we would have change from type-II SC to type-I SC 
near TmQ. This result is apparently incorrect in the RE ternary compounds, 
since 6ex ~ Gem, thus making p finite and k stays practically unchanged. So, the 
change of the type of transition near T^q is not due to the change of k but it is 
due to the much faster temperature falloff of Hc2{T) than of H(.{T). We shall 
not discuss further this interesting subject but refer the reader to [3] where 
various phases in the H-T phase diagram were analyzed. Depending on the 
demagnetization effects several phases can be realized in the same material, 
such as Meissner-, vortex-, LOFF- or even intermediate-phase. 



5 Josephson effect on magnetic superconductors 



After the remarkable theoretical discovery by Buzdin and coworkers of the 
possibility of tt- Josephson junctions in the hybrid S — F — S structure where 



20 




Figure 3. The Josephson junction with the insulating contact. Sl and Sr are su- 
perconductors with spiral magnetic order. The exchange fields hi^R at the surface 
make angles ^^^/jwith the y-axis. Ql,r are along the z-axis. 

F is a ferromagnet [9], [10] the interest in Josephson junctions with magnetic 
degrees of freedom has grown dramatically - see this issue. In that sense it is 
natural challenge to investigate this problem in magnetic superconductors. 



5.1 -n -contact due to triplet amplitude F^^ (Fj 



[1/ 



In [7] the tunnelling Josephson junction was studied with the left-L and right- 
R bulk magnetic superconductors in which the spiral magnetic ordering is 
realized - see Fig.3. The spiral magnetic order is characterized by the wave 
vector Ql,r and the exchange fields hg^^^^^ = hL{R)e^^'^'^^\ respectively, while 
superconductivity in the banks is described by the order parameter A^^r =| 
^L,R I e**^^'-^. Due to simplicity it was assumed | | = | A/j |= A, hi = hR = 
I Ql 1 = 1 Q_R 1= Q where Ql,b. = Xl,rQ^ are orthogonal to the tunnelling 
barrier with the spiral helicity X-L(R) = il for Ql,r along (+) or opposite(-) to 
the z-axis. Note, that such a junction is characterized by the superconducting 
phase (p = (pL — '^R 7^ and magnetic phase 9 = 9l — 9r ^ 0. It turns 
out that besides the singlet amplitude the triplet pairing amplitudes 

(Fii) play very important role in the Josephson effect [7]. (Note, that 
-F|| (F^) was first introduced in [4] in the study of the spiral magnetic order 
in SC and first applied to the Josephson junction in [7]. Later on this effect 
was rediscovered by the Efetov's group in studying S-F-F-S structures with 
rotating magnetization [26], where (F^) give rise for new effects see this 
issue.) It turns out that the Josephson current contains two parts 



21 



Jj{^, 0) = {J' - Ji^ cos e) sin (p, 



(55) 



where 

r^r^T Y: \T^„^,\'F{^{kL,u;n)Fn{kR,-u;n) (56) 
(and ~ A^) is due to the singlet amplitude and 

kz,,ki{,a;„ 

+ Fl{kL,uJn)[Fl{kR, -ujn)]*} (57) 

is due to the triplet amplitude. It turns out that 

Ji^ ~ A^/i^i/, + (XLX«)/2(A, /i)], (58) 

where /i^2(A, /i) are given in [7], while x = XlXr is the total helicity (in [26] 
renamed to chirality) of the junction. It was shown in [7] that in some param- 
eter region the triplet effects dominate, i.e. | Jt^ \> J*', thus giving rise to the 
TT-Josephson junction. Prom Eq.(54) it is clear that by changing the magnetic 
phase 9 and chirality x one can tune the system from 0— to tt— junction. This 
new degree of freedom in the junction - the magnetic phase 6, first proposed 
in [7], opens a new possibility for switching elements and quantum computing. 
Prom the physical point of view the above model is a paradigm for analogous 
effects in S-P-P-S structures with rotating magnetization, in which case is 9 
the angle between magnetizations in neighbouring layers. 



5.2 Combined superconducting and magnetic Josephson effect 



In [8] the above model is developed further by including the tunnelling of 
electronic spins and their effect on the energy of the contact. Namely, in fer- 
romagnetic superconductors with rotating of magnetization (such as spiral) 
besides the standard Green's function G|| (G^), F|j^(F||) for singlet SC other 
Green's functions G-^i (G^j,t) ^^"^ -^TT (-^iJ.) important [7], [8], since they 
can produce a static spin current Jspin through the junction (in absence of 
voltage) , 

Jspin J spin, G ~l~ J spin, F (^^) 
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where 

Jsp^n,G -Y.\T? {G^i,LG^^,R - G^^,lG^^,r] ^ h'sinO. (60) 

Jspin,F ~ 5^ I T |2 [F^^{kL,UJn)[F^^OiR, -UJn)]* 

- Fl{kL,uJn)[Fl{kR,-uJn)Y} h'A^f^cosipsine. (61) 

The exact expression for Jgpin.G and Jspin,F wiU be pubhshed elsewhere [8]. 
The energy of this combined magnetic and superconducting Josephson junction 

Emj{e) + Ej{ip,e) is 

E{9,ip)^ -Ah^ cos ^ - (S + C^h^ cos 6) cos ip, (62) 

where the cumbersome expressions for A, B,C are given in [8]. Note, that both 
the spin Jspin{0, ~ dE/dO and the charge J,/(v9, 9) ~ dE/dip Josephson cur- 
rent depend on tp and 6. Thus, by tuning 6 and p one can tune these currents. 
In case of small contacts with small charge and "spin" capacitance the system 
is in the quantum regime thus giving possibility for a novel Josephson qubit. 
In fact the latter consists from two qubits - the superconducting and magnetic 
one, which is of a potential interest for applications [8] . 



6 Conclusion 

The rare earth ternary compounds arc reach physical systems which allow 
coexistence of superconductivity and various magnetic orders, such as ferro- 
magnetic, antiferromagnetic, weak ferromagnetism. It turns out that in these 
systems superconductivity and ferromagnetism never coexist and the latter is 
modified into a spiral or domain structure - depending on magnetic anisotropy. 
This is realized in rare earths ternary compounds as well as in Auln2 where the 
modified nuclear ferromagnetism and superconductivity coexist. Although the 
antiferromagnetic order and superconductivity coexist much easier, these sys- 
tems show peculiar behavior in the presence of nonmagnetic impurities which 
surprisingly act as pair-breakers. In case when the antiferromagnetic order is 
accompanied by the weak ferromagnetism new coexistence phases appear - 
the Meissner and spontaneous vortex state. Magnetic superconductors show 
peculiar behavior in magnetic field. Near the magnetic critical temperature 
the upper critical field goes to zero faster than the thermodynamical field, 
thus giving rise to the first order transition. Various phases are possible in the 
H — T diagram depending on the purity and demagnetization effects of real 
samples. The lower critical field is weakly affected by the exchange field due 
to localized moments. 
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The Josephson junctions based on bulk ferromagnetic superconductors with 
spiral order are characterized by the superconducting and magnetic phase, 
opening possibilities for a new kind of coupled qubits. The triplet pairing 
amplitude gives rise to the tt— junction which can be tuned by changing the 
magnetic phase and chirality. 

This article is submitted for the Special Issue Comptes de I'Academie des 
Sciences: Problems of the Coexistence of Magnetism and Superconductivity 
edited by A. Buzdin. 
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7 Appendix 

The problem of the coexistence of SC and magnetic order with a wave vector 
Q can be studied in principle by using Gorkov equations for any Q. However, 
in systems where a « is fulfilled the quasiclassical Eilenberger-Larkin- 
Ovchinikov (ELO) equations are more suitable and efficient. 

7. 1 Gorkov equations for MS 

These equations contain normal and anomalous Green's functions {G)ap{x, y) = 
- < Ti'a{x)'i/jl{y) > and {F^)ap{x,y) =< f'ijjl{x)'iljl{y) >, where x = if.r) 
and a, (5 =t, [. The superconducting order parameter is defined by 

A*(r-; - = -\v{r, - f,)[< >-< >\ (63) 

where the pairing interaction V{fi—rj) is responsible for the s—wave pairing in 
absence of magnetic order. Note, that in the presence of any (in) commensurate 
magnetic order the superconducting order parameter should be nonuniform, 
i.e. A(ri,r2) = ^(^i ~ ^2,'^, ^ = {'^i + o)/^' However, it has been shown 
in [4] that when h < hvpQ the nonuniform part of A is small and of the order 
{h/hvpQ)^- That is the reason for our assumption A*(rj,rj) = A*(rj — f^). 
The set of equations for the Green's functions is given by 

- eo(p) - Ve.{^\G.Sr. ^ + / ^^'^A(f - f)F^Jf , r') = 5{r - r^). 
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K + e^o(p) + f") - I d'xA*{f- x)G^„{x, r-) = 0, (64) 



where 



A = A 


1 










-1 





(65) 



In case of a spiral magnetic ordering (with spatial rotation of magnetization) 
besides the singlet pairing amplitude ~ A there is also a triplet pairing 
amplitude F^i^ ~ A • /ig (and F^) [4]. is responsible for singlet pairing 
with the order parameter A|j(r), while the triplet amplitude which is 
due to the rotating magnetization which mixes spin up and down, does not 
give rise to the triplet pairing since it is assumed from the very beginning that 



HI 



0. 



Note that Fj^ 



contains two terms - the first one is odd in 



frequency cun and the second one is odd in momentum k. 



7.2 ELO equations for MS 



In problems related to the presence of the exchange field h(r) = h{r)ez (ori- 
ented along the 2;-axis) and of nonmagnetic potential and spin-orbit scattering 
a generalization of the ELO equations is needed [12], [27]. In that case the 
Gor'kov equations contain 4x4 Green's functions G{xi, X2) = — (T^(a;i)^^(a;2)) 
with the four-component spinor ^^(xi) = ('0|'0|'0tV'i)- Here fj are Nambu ma- 
trices in the electron-hole space and o"j are Pauli matrices in the spin-space 
and X = (r. r). (Note, that here the functions F are defined with minus sign.) 
The SC order parameter is A(r) = if_|_(72A(r) — if_o"2A*(r) with 

A(r) = iVTr{r+a2G(x, x)}/2. (66) 

By integrating out unimportant degrees of freedom at small distances one 
defines the quasiclassical Green's function g^^CR^pp) 

gu,M:PF) -jd^pj d{ri - r2)e-^P«"-i-'-^)f3G'(ri,r2,a;n). (67) 

Then the ELO equations for the quasiclassical Green's functions read 

iypVngo.^ = [tujr,h + fsA(R) + h(R)as - f3EX^(R, p^), ^.J, (68) 
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where the impurity self-energy in the Born approximation is given by 

KZ'i^^Pp) = c,N{0){U{pj. - p'j.)9^AR,p'f)U{pf - p^))p^. (69) 



The impurity potential (matrix) contains the non-magnetic and spin-orbit 
scattering 



U{Pf - p'f) = Uih + iUso[pF X p'f\&, 



(70) 



with the vector matrix a = [(1 + fs)a+ (1 — f3)cr2cr(72]/2 and = pp/ \ pp 
The matrix g^^^ (R, p^) is given by 



^9+ iU^ 
g- if- o 
-ifl -g. 
[-if I -g+) 



(71) 



By including also orbital effects in magnetic field, which is determined by the 
vector potential A(R), the ELO equations written in components g± and f± 
read 



pn,± + ievp ■ A(R) ± ih{R) - -Vp ■ Vr]/±(pf, R, u^n) 
= A±(pF,R,u;n)g'±(pir,R, a;„), 

[Cjn,± + ie^F ■ A(R) ± i/i(R) + ^vj. • Vr]/1(pf, R, w„) 
= A^(pir, R,a;„)5(±(pF, R, a;„), 



(72) 



(73) 



with the normalization condition and the self-consistency equation, respec- 
tively 



gl + fif^ = l 



(74) 



A(p^,R) = — V^(PF,P^){/+(R,PF,'^n)+/-(R,P^,^n)}. (75) 



Cjn,± and A±(pi;', R, a;„) are defined by 
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+ :^{g^{p'F, R, a;,) sm'{9 - ^'))p^, (76) 

r 

A±(pF,R,a;n) = A(pj.,R) + y(/±(p^,R,a;„))p^ 

+ ^soiUiPF, R, ^n) sin2(^ - ^')>P^ (77) 



where Ti = CiN{0) \ Ui \^ and = QiV(O) | 



|2 



The microscopic theory of magnetic superconductors which takes into account 
spin (exchange) and orbital (electromagnetic) effects of magnetic order, as well 
as of nonmagnetic impurity scattering, has been developed by using the above 
generahzed ELO equations [3], [12]. 
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